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ABSTRACT 


The application of the 2D-video disdrometer to measure fall speed and snow size distribution 
and to derive liquid equivalent snow rate, mean density-size and reflectivity-snow rate power 
law is described. Inversion of the methodology proposed by Bohm provides the pathway to use 
measured fall speed, area ratio and ‘3D’ size measurement to estimate the mass of each 
particle. Four snow cases from the Light Precipitation Validation Experiment are analyzed with 
supporting data from other instruments such as Precipitation Occurrence Sensor System 
(POSS), Snow Video Imager (SVI), a network of seven snow gauges and three scanning C- 
band radars. The radar-based snow accumulations using the 2DVD-derived Z e -SR relation are 
in good agreement with a network of seven snow gauges and outperform the accumulations 
derived from a climatological Z e -SR relation used by the Finnish Meteorological Institute (FMI). 
The normalized bias between radar-derived and gauge accumulation is reduced from 96% 
when using the fixed FMI relation to 28% when using the Z e -SR relations based on 2DVD data. 
The normalized standard error is also reduced significantly from 66% to 31%. For two of the 
days with widely different coefficients of the Z e -SR power law, the reflectivity structure showed 
significant differences in spatial variability. Liquid water path estimates from radiometric data 
also showed significant differences between the two cases. Examination of SVI particle images 
at the measurement site corroborated these differences in terms of unrimed versus rimed snow 
particles. The findings reported herein support the application of Bohm’s methodology for 
deriving the mean density-size and Z e -SR power laws using data from 2D-video disdrometer. 
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1. Introduction 


The measurement of liquid equivalent snow rate (SR) from radar has long been 
recognized as a difficult problem in quantitative precipitation estimation (QPE) but one 
of great importance for weather forecasting, hydrology, detection of aviation hazards 
and other remote sensing applications (e.g., ground validation for microwave radiometry 
from space). The validation of QPE by radar is difficult at best given the fact that 
accurate measurement of winter precipitation by gauges remains challenging due to the 
sheer variety and variability of physical properties which can change dramatically with, 
for example, relatively small changes in environmental conditions. Some of the 
important physical properties that one could list, for example, are (i) ‘3D’-size, (ii) 
terminal fall speed, (iii) particle size distribution, (iv) density (or, mass), (v) shape, (v) 
composition and (vi) porosity. Some of these attributes are not independent as 
evidenced by the large literature that exists in describing density (or, mass)-size and fall 
speed-size relations for different kinds of winter precipitation (e.g., Pruppacher and Klett 
2010; Mason 2010). The fall speed is also dependent on shape, composition and 
porosity. Thus, it follows that fall speed is fundamental to characterization of frozen 
precipitation followed by a good measure of ‘3D’-size, particle size distribution and 
porosity. From the radar reflectivity perspective, the ‘3D’-size and associated size 
distribution and the density (or, mass)-size relation is of fundamental importance. For 
Rayleigh scattering the reflectivity is directly related to E [m(D) 2 ] where m is the particle 
mass and E stand for expectation or integration over the size distribution (Ryzhkov et al. 
1998); however, the mass is not easily measured on a particle-by-particle basis. On the 
other hand, the liquid equivalent snow rate (SR) is directly related to E [m(D) Vf(D )] 
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where Vf is the fall speed. It follows that empirical Z e -SR power laws can be derived if 
mass-D and Vf-D power laws are assumed (e.g., Matrosov et al. 2009) and the size 
distribution is measured (e.g., Sekhon and Srivastava 1970). The simulations of 
Matrosov et al. (2009) suggest that the overall uncertainty of estimating SR from 
reflectivity measurements can be as high as a factor of 3 or so. A more direct method is 
to correlate Z e from radar with SR measured by snow gauges (e.g., Fujiyoshi et al. 1990 
and references therein) which can lead to climatological Z e -SR power law. The advent 
of optical-based surface disdrometers, however, has led to more accurate methods to 
characterize the physical properties of snow, leading to m-D, Vf-D and area ratio-D 
relations that are consistent via hydrodynamic theory (Bohm 1989; Mitchell 1996; 
Heymsfield and Westbrook 2010). Combined with scattering models (size, shape, 
dielectric constant), it leads to more consistent Z e -SR power laws (Huang et al. 201 1). 


There are a number of disdrometers (mainly optical) that are available (some 
commercial and others in the research category) that measure a sub-set of the physical 
parameters listed above (only instruments that can image the particles are considered 
here). Hanesch (1999) and Schonhuber et al. (2000) used the 1 st generation 2D-video 
disdrometer (2DVD) which measures fall speed and two orthogonal images from which 
an apparent volume (also, size distribution based on ‘3D’-size) as well as an estimate of 
porosity (via the area ratio to be described later) can be computed. Later, Brandes et al. 
(2007) used the 2DVD to estimate the coefficient and exponent of a mean density-Do 
power law (mainly for fluffy snow aggregates) by comparing 15-min liquid water 
accumulations with a collocated Geonor gauge (Do is the median volume diameter of 
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the particle size distribution). They also examined the particle size distribution in detail 
by fitting with a gamma model and deriving correlations between the model parameters 
(e.g., shape parameter ;u and slope parameter A). They conclude (in their Section 6, 
page 648) that, “...The video disdrometer is a powerful observational tool for studying 
the microphysical properties of winter storms”. Further, Brandes et al. (2008) 
investigated power law relation between terminal fall speed and size and its 
dependence on temperature. The use of radar and 2DVD for estimating density-size 
and Ze-SR power laws is described by Huang et al. (2010; 201 1) whereas Zhang et al. 
(2011) demonstrated the importance of density-size power laws (empirically adjusted by 
fall speed) in comparing 2DVD-based reflectivity with ground radar. While the 2DVD is 
commercially available, a similar research instrument HVSD (Hydrometeor Velocity Size 
Detector; Barthazy et al. 2004) measures the fall speed and projected image in one 
plane. It has been used by Zawadzki et al. (2010) to investigate the natural variability of 
snow terminal velocity with size. They concluded that the exponent of the terminal 
velocity-D power law could be fixed at 0.18, while the coefficient is variable from event- 
to-event. Szyrmer and Zawadzki (2010) describe a methodology to derive the average 
relationship between terminal fall velocity and the mass of snowflakes via elaboration of 
the methodology of Bohm (1989) proposed by earlier Hanesch (1999); the latter used 
the 1 st generation tall 2DVD design. In fact, the development of the HVSD by Barthazy 
et al. (2004) followed Hanesch and lead to a simpler instrument with two parallel light 
planes but with much slower line scan frequency camera. The work described herein 
follows Szyrmer and Zawadzki (2010) but uses the 2 nd generation low profile 2DVD 
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(Schoenhuber et al. 2008) to derive the Z e -SR power law with validation provided by a 
network of seven snow gauges. 


Another research instrument is the Snow Video Imager (SVI; Newman et al. 2009) 
which, unlike the line scan camera, uses a CCD (charge-coupled device) full frame 
camera (60 frames per second) and images are obtained almost simultaneously; 
however, it does not measure the fall speed. SVI software yields a size estimate of each 
particle as an equivalent diameter that corresponds to a circular equivalent-area 
diameter of the irregular shape (with holes filled). Some advantages of the SVI over the 
2DVD is that it has a large sample volume (twice that of the 2DVD), better pixel 
resolution (nominally 0.05 mm by 0.1 mm) and its measurements are less sensitive to 
wind. In this work, the SVI is mainly used to determine the particle size distribution for 
comparison with the 2DVD, and to examine samples of images to distinguish between 
unrimed and rimed snow particles. A new commercially available instrument is the Multi- 
Angle Snowflake Camera (MASC; Garrett et al. 2012) which gives high resolution (IQ- 
50 pm) photographs of snow particles from three viewing angles, along with their fall 
speed. One disadvantage is that the sample volume is small (about 1/10 of the 2DVD). 


This article is organized as follows. In Section 2 the specific details of estimating the 
apparent volume and ‘3D’-apparent diameter, the adjusted particle size distribution and 
the application of Bohm’s (1989) method are described. Section 3 constitutes the main 
bulk of the article and describes the 2DVD processing and derived products culminating 
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in Ze-SR power laws for the four snow days, comparison of liquid equivalent snow 
accumulations derived from radar with a network of 7 snow gauges, and radar-based 
accumulation maps. A short summary and conclusions are given in Section 4. 

2. The basis for snow measurements using the 2D-video disdrometer 

2. 1 The apparent volume and diameter 

The 2DVD gives two views (front and side views; actually silhouettes) of the particle in 
two orthogonal planes as shown in the example in Fig. 1. It is obvious that the ‘true’ 
volume of such an irregular particle cannot be calculated and thus we define here the 
apparent volume ( VL app ) assuming that the particle is an ellipsoid. The apparent volume 
is defined as an average of two ellipsoidal volumes: 

Vl app = - b D app * = VLa ^ La vv 2 (1) 

where the apparent diameter is D app and, 

VL appl =l(H*Wi*W 2 ) (2) 

where, 

H = ^H 1 * H 2 (3a) 


Wi,2 = 


4*A e i,2 

n*H 


(3b) 


The A e i ,2 are the shadow areas (see Fig. 1) from the two views. The second ellipsoid 
estimate, VL apP 2 , is defined as: 

VL app 2 W maxl * W max2 ) (4) 
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where, 


HH = JHHj. * HH 2 (5a) 


HHi2 = ^±^ 

' Tl*W max i2 


(5b) 


In (5b) the Wmax equals the maximum width of the scan line or ‘slice’ (measured from left 
to right in Fig. 1); the subscripts 1,2 refer to maximum width as determined from each 
view. The method of calculating VL apP and D app here generally follows Hanesch (1999) 
which is somewhat different from Schonhuber et al. (2000) which was used later by 
Brandes et al. (2007). Also, the apparent diameter ( D app ) is different from the ‘size’ 
measured by instruments that give the particle image in only one plane such as aircraft- 
mounted imaging probes (which give the top view). The ‘size’ is often defined as the 
maximum distance between two pixels or the diameter of the smallest circle that 
completely circumscribes the image or the equivalent-area diameter (Hogan et al. 
2012). The latter also define the mean diameter as the mean of the particle dimensions 
in two orthogonal directions which they found to be better related to radar reflectivity. 
Since the true volume of snowflake is not known, the accuracy of our method of 
calculating VL app cannot be determined. However, from the simulations of Wood et al. 
(2012) who used ellipsoidal shape models with canting it is can be inferred that the 
apparent diameter defined here gives a more ‘realistic’ measure of ‘3D’ size made 
possible by the availability of two orthogonal images from the 2DVD. 


2.2 Snow size distribution (SSD) 
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In a certain time window (typically 60 seconds for 1 -minute averaged size distributions), 
all ‘matched’ snow particles are sorted into M size bins according to the apparent 
diameter ( D app ) and the ‘un-adjusted’ size distribution N m (Di) is computed as: 


NmW 


1 yiA/^i 1 

At*A£> "/= 1 Aj*vj 


[mm 1 m 3 ] (6) 


where D, is the center diameter of the i th size bin (from 1 to M) in mm; AD is the bin 
width in mm; A\ is the measurement area in mm 2 ; Vj is the fall speed in m s' 1 and At is 
the time window in seconds. The fall speed measurement is fundamental to the 2DVD 
and relies on the ability to match the particle that falls in the upper light plane (and is 
imaged by Camera A) to the same particle that falls through the lower light plane and is 
imaged by Camera B (see Fig. 2). The match criteria used here are adapted from 
Hanesch (1999) as elaborated by Huang et al. (2010). If the match criteria are not 
satisfied then that particle is rejected; it follows that the concentration will tend to be 
under-estimated. To re-adjust the measured N m (Dj) for this underestimate (assumed to 
be a constant factor /) the following procedure is used. 


Assume that snow falls uniformly over the instrument. Then, the theoretical number of 
snowflakes falling through the virtual measuring area divided by the theoretical number 
of snowflakes falling in the scan area of each camera (shown in Fig. 2) should be equal 
to the ratio of these two areas as: 


theoretical # of snowflakes in virtual measuring area 
theoretical # of snowflakes in scan area of single camera 


100 « 
= 0.4 


250 


( 7 ) 


Therefore, an adjustment factor ^is derived as: 
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( 8 ) 


_ 0.4*(# of snowflakes actually counted, in scan area of single camera ) 

^ # of mached snowflakes in virtual measuring area 

The “re-adjusted” concentration in each size channel (N(Di)) is defined as: 

JV(A)=y*J\UD0 (9) 

where yis assumed constant (y> 1). In essence, the “raw” or unadjusted SSD is simply 
scaled by the factor y The validity of this adjustment will be evaluated by comparison 
with SSD from the snow video imager (Newman et al. 2009) as well as determining the 
y independently by comparison with the SVI as described later in Section 3.1 . 

2.3 Bohm’s Method 

Bohm (1989) developed a general methodology for the terminal fall speed of solid 
hydrometeors based on the mass, the mean effective projected area ( A e \ see Fig. 3, 
also referred to as shadow area) presented to the flow, and the smallest circumscribed 
area (A\ circle or ellipse depending on the shape of the snow particle). Since the 2DVD 
can measure the fall speed of each snowflake as well as two orthogonal images (Fig. 
1), we are able to compute the mass of each snowflake by inverting the Bohm 
equations. The assumption is that A e which is the projected area in a plane normal to 
the flow, is approximately equal to the area from the side or front views (measured by 
the 2DVD) for irregular shaped particles. This assumption has been evaluated as being 
‘reasonable’ by Szyrmer and Zawadzki (2010) who use the HVSD which gives the side 
view only. We first compute the Reynolds number (Re) from fall speed (Vi) and viscosity 
(rj) as: 
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Re 


2 * p a * V f /Ay 
rj \nJ 


(10a) 


where the characteristic dimension is the ‘area’ diameter (A being the area of the 
smallest circumscribed ellipse or circle that completely encloses the particle image). 
The r/ and air density (p a ) are computed from temperature, air pressure and humidity. 
Next, we compute the Davies number (X) from Re as: 


X = 


kt)^ 

2 ) 
- 1 

0.1519 

V 

> 


(10b) 


Finally the mass of the snowflake is computed as: 


n *rj 2 *X /7l e \ 1/4 
8*g*p a \AJ 


(10c) 


where g is the acceleration due to gravity. The ratio ( AJA ) is referred to as the ‘area 
ratio’ or A r which is < 1. The MKS units are appropriate for the variables in (10). The 
relative error in the estimate of mass due to uncertainty in the fixed relation between X 
and Re, and in the estimation of A r has been evaluated by Szyrmer and Zawadzki 
(2010) as between 40-50%. The propagation of error from (10a) to (10c) is complicated 
and the reader is referred to the aforementioned reference for details. 


The calculation of the minimum circumscribed area (A) is based on the rectangle which 
completely encloses the particle (the rectangle width is 1/14 and height is H; see Fig. 1). 
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We first assume that, (i) A is the maximum ellipse that can be fitted inside the rectangle 
and compute the area ratio A e /A which should be < 1 . If this ratio is greater than 1 , we 
assume that, (ii) A is the minimum circle that can contain the rectangle. The minimum 
circumscribed area estimated from (i) usually tends to underestimate A whereas from 
(ii) tends to overestimate A. The apparent volume ( VL app ) and apparent diameter ( D app ) 
were defined earlier, thus the density ( p ) is obtained as the ratio of mass (m) to VL app for 
each particle. Since our measurements are restricted to frozen ice precipitation, the 
density is also restricted to rr\\r\[mA/L app 0.917] in cgs units. The mean density is 
calculated for each size bin and a power law fit of the form p=aD app 13 is obtained for the 
precipitation event. Here D app is in mm and p is in g cm -3 . The mass-D a pp power law then 
is m=ct(n/6)D app P +3 . 

For an area sampling measurement device such as the 2DVD, the liquid equivalent 
snow rate ( SR m ) can be computed directly as: 


where N is the number of size bins, M is the number of snowflakes in the i th size bin in 
At (typically 60 seconds), VLj is the liquid equivalent volume of j th snowflake in mm 3 (this 
is the product pVL app where p as a function of D app is given by the power law fit) and Aj 
is the measurement area for the j th snowflake in unit of mm 2 . The adjusted snow rate is 
SR- y *SR m where the /-factor was defined in Section 2.2. 


N M 



i= 1 j = 1 


[mm hr x ] (11) 
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We use the T-matrix (Waterman 1971; Barber and Yeh 1975) method to compute the 
radar cross-section of each particle and the equivalent reflectivity assuming: 

• Refractive index: computed by the Maxwell-Garnet (1904) mixing formula 
with temperature from environmental data. The particle is assumed to be a 
mixture of ice inclusions within an air matrix with effective density p as a 
function of D app as given by the 2DVD-derived power law fit 

• Orientation: the zenith angle is Gaussian with zero mean and 45° standard 
deviation; the azimuthal angle is uniformly distributed in [0 tt] 

• Particle Size Distribution: as in Section 2.2 which defines N(D)=y*N m (D) 

• Particle shape: oblate spheroid with volume = VL app and axis ratio = 0.8 

3. Experimental Data from LPVEx 

The Light Precipitation Validation Experiment (LPVEx) was held in the Fall of 2010 in 
the area surrounding Helsinki, Finland as a collaborative project between the Finnish 
Meteorological Institute (FMI), University of Helsinki and the NASA Ground Validation 
program (Petersen et al. 2011). While the experiment had multiple objectives and 
extensive suite of instruments, the focus herein is on snow measurements made at the 
Jarvenpaa site with the 2D-video disdrometer; this site also had the Precipitation 
Occurrence Sensor System (POSS; Sheppard and Joe 2008), the Snow Video Imager 
(SVI; Newman et al. 2009) and an OTT-PLUVI02 gauge with Tretyakov and Alter 
shields (lanza et al. 2006; Rasmussen et al. 2012). Three C-band polarimetric radars 
located at Kumpula, Vantaa and Kerava (Koskinen et al., 2011) provided for excellent 
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coverage over the Jarvenpaa site as well as over the network of six FMI snow gauges. 
Fig. 4 shows the location of the 3 radars, the Jarvenpaa site and the gauge network. 
Briefly, the three radars are nearly identical with 1° beams and using simultaneous H-V 
polarization on transmit and simultaneous reception of the H and V polarized back- 
scattered signal components via two receivers. The minimum detectable Z e is about -10 
dBZ at range of 50 km. The reflectivity data from each radar covering matched areas of 
precipitation were used to construct the CDF of Z e enabling accurate cross-calibration of 
the radars (Hirsikko et al, 2013). All the radars are Vaisala dual-polarization weather 
radars, a detailed description of Vantaa radar operations is presented by Saltikoff and 
Nevvonen, (2011). 


Table 1 lists the four snow days where there was significant precipitation in Helsinki and 
surrounding areas. The snow events on these days were also favorable for 2DVD and 
other snow measuring instruments as the wind speeds were < 4 m s' 1 at the Jarvenpaa 
site. As seen in Table 1, the 30 Dec 2010 case could be sub-divided into two snow 
events based on the synoptic conditions. Similarly, the first event on 12 Jan 201 1 (0800- 
1230 UTC) could be separated from the second event that covered the period 2230- 
2359 UTC which further continued the next day (13 Jan) until 0500. The liquid 
equivalent snow accumulations ( SA ) from the OTT-PLUVIO gauge ranged from 1.5 to 
4.2 mm. 


3 . 1 Example of 2D VD processed data from 30 Dec 2010 
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As mentioned earlier, one of the fundamental measurements provided by the 2DVD is 
the fall speed, an example of which is provided in Fig. 5 from the first event on 30 Dec 
2010. The instrumental error in measuring terminal fall speeds is < 4% (for fall speeds 
<10 m s' 1 ; Schonhuber et al. 2008). Such high accuracy is due in part to the plane 
distance calibration which is performed frequently and accounts for slight deviations in 
the plane distance depending on the location within the virtual measurement area (see 
Fig. 2); further the line scan frequency is quite high close to 55 kHz. Zawadzki et al. 
(2010) evaluated the fall speed measurement error for the HVSD which, to the best of 
our knowledge, does not account for plane distance deviations within the measurement 
area plus the line scan frequency is much lower, closer to 10 kHz. They estimated that 
the instrumental uncertainty for the HVSD is around 12% for fall speeds below 2 ms' 1 . 
While a similar analysis has not been done for the 2DVD, the contribution of 
instrumental error to the fall speed is expected to be much smaller than the natural 
variability which is depicted by the ±1o bars in Fig. 5. It is also evident that the 
commonly used power law fit for Vf versus D, while analytically convenient, does not fit 
the data as well as an exponential fit of the form \4= c[1 - d*exp(- KD app )]. 


The snow size distribution (SSD) for the same event is shown in Fig. 6a where the 
distribution from 2DVD is compared with that derived from the Snow Video Imager 
(SVI). The -factor was estimated as 2.21 (Section 2b and eq. 8). Fig. 6b shows similar 
plot for 6 Jan. 201 1 event. The agreement in the SSD is quite good given that the two 
instruments are based on distinctly different measurement principles and sample 
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As a further check on the estimation of y using (8), the unadjusted SSD from the 2DVD 
has been forced to match the SVI in each size bin and the resulting mean ysvi is 
computed as: 


where N is the number of size bins, Nsvi(Di) is the SVI-measured concentration for the 
i th bin, and the corresponding 2DVD-measured N m (Di) is obtained as in (6). For the case 
shown in Fig. 3.3 the ysvi was found to be 2.46 which is in close agreement with y= 
2.21 . For the other snow events listed in Table 1 the /comparisons are given in Table 4. 
As noted in the introduction the SVI gives a measure of the equal-area circular diameter 
which is not the same as D app from the 2DVD. We ignore the different estimates of D 
from the two instruments is so far as validation of the single camera-2DVD based /- 
factor estimation is concerned. A more elaborate discussion of SVI estimation of 
different measures of D and related characterization of uncertainties in estimation of Z e 
and SR are given in Wood et al. (2013). 

The area ratio ( A r ) discussed in Section 2c plays an important role in inverting Bohm’s 
methodology to derive mass from the fall speed. Schmitt and Heymsfield (2010) 
comment that, “...area-dimensional and mass-dimensional relationships are rarely 
developed from the same dataset”. Fig. 7 shows the frequency of occurrence plot (in log 
scale) of Ar vs. Dapp for the same 30 Dec 2010 snow event. Also shown are the bin 
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( 12 ) 


averaged mean and ±1a standard deviation bars along with the power law fit 
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A r = 0.71 Dapp' 008 . The variability in A r is quite large but in general agreement with 
Zwadzki et al. (2010) who used data from the HVSD but allowed A r > 1 . The mean fit in 
Fig. 7 is in good agreement with that given in Zawadzki et al; they obtain A r =0.75D~° 17 
(but their ‘D’ is the maximum dimension from the side-view image). A somewhat 
different power law fit was obtained by Schmitt and Heymsfield (2010), who used cloud 
imaging probe on aircraft penetrations of ice clouds aloft (this is not surprising since our 
results are at the surface in heavier snowfall). Schmitt and Heymsfield obtained an 
exponent of -0.25 for the ARM data set (Heymsfield et al. 2004), but their coefficient 
was lower by a factor of 2. 

The final result from Bohm’s methodology is the ability to derive a mean density-D app 
power law and Fig. 8 shows the same for the 30 Dec 2010 event. While there is large 
variability in density for a given D app (especially evident for small particles D app < 1 mm 
which might be related to difficulty in matching such particles from the two camera 
images resulting in erroneous fall speed determination); nevertheless, there is an 
inverse relation between density and D app and the power law fit is p=0.15 D app °- 86 for 
this event (Table 1 gives the coefficient and exponent for the other events). Plots of p 
versus D app from Table 1 using the coefficient (a) and exponent () 3 ) found herein for the 
four snow days are close to the mean climatological relation found by Brandes et al. 
(2007; a=0.178, /3=- 0.922) as well as Holroyd (1971; cr=0.17, /3=-1 ) and Fabry and 
Szyrmer (1999; o=0.15, /3=-1) with the caveat that ‘D’ in each of the quoted references 
are not calculated in the same manner (see, also, Table 2 from Brandes et al. 2007). 
The exponent of the mass-D app power law is given by 3+/3; from Table 1, the latter 
exponent varies between 2.04 to 2.21 generally within the range obtained by Schmitt 


342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 


and Heymsfield (2010) based on fractal simulations of large aggregates (range between 
2. 1-2.2), and close to the experimentally obtained exponent of 2.2 for the ARM dataset 
(Heymsfield et al. 2004). 


The liquid equivalent snow rate (SR) for the 30 Dec 2010 is calculated as given in (11) 
using the mean p-D app power law fit from Table 1 for the two snow events that occurred 
on that day. The SR is adjusted by the / factor. Fig. 9 shows the liquid equivalent snow 
accumulation from the 2DVD compared with the collocated OTT-PLUVI02 gauge at the 
Jarvenpaa site. The maximum SR during the two snow periods occur at around 1230 
and 2200 UTC. The agreement between 2DVD and gauge is very good for this event 
(accumulations are based on 1-min SR from 2DVD). From Table 1, the accumulations 
between the 2DVD and gauge for the other days are also in good agreement. It is 
difficult to estimate the accuracy of the 2DVD-derived snow rate but assuming the p- 
Dapp is ‘exact’ and valid for the entire event, the dominant systematic error would be in 
the ^adjustment parameter of the SSD. Otherwise, systematic error would primarily 
arise due to incorrect estimate of a and secondarily /3. 

3.2 Reflectivity and Z e -SR power law 

The reflectivity ( Z e ) at C-band (frequency 5.5 GHz) is computed from 2DVD-measured 
1-min averaged N(D app ) and the mean p-D app power law fit (for the entire event), based 
on the assumptions listed towards the end of Section 2.3. It is well established that for 
Rayleigh scattering and using the Maxwell-Garnet mixing formula (ice inclusions inside 
an air matrix) that Z e can be expressed as: 
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/ 1 \ 2 I I 2 r D max 

Z e = PsnowD 6 N(D)dD ( 12 ) 

where |K /ce ,w| 2 are the dielectric factors of solid ice and water. Since the mass of the 
particle is m=pVL a p P , it follows that Z e can be simply computed (suppressing constants) 
as the sum(/77 2 ) over all the particles. Thus, the reflectivity is very sensitive to the m-D app 
relation (or, equivalently the p-D app ) as shown by a number of previous studies (e.g., 
Ryzhkov et al. 1998; Matrosov 2009). Errors can arise from uncertainty in the ^factor 
which scales the N(D) or uncertainty in the coefficient a and less so in the exponent /?. 
Note that the T-matrix scattering code is used to compute Z e at C-band frequency. 

Fig. 10 shows time series comparison of Z e from 2DVD, POSS [at Jarvenpaa site for (a) 
30 Dec 2010 event and (b) 6 Jan. 2011] and the scanning Kumpula C-band radar 
reflectivity data extracted over the same site (areal average over 1°X 1 km). The 2DVD 
and POSS reflectivities are 1-min averaged whereas the Kumpula radar data were 
available every 5 min. The 2DVD data are somewhat more ‘noisy’ as compared to 
POSS due mainly to sampling error (the POSS has a very large sample volume by 
several orders of magnitude relative to the 2DVD). The sampling error in the 2DVD 
measure of Z e was evaluated by Huang et al. (2011) by using two 2DVD units located 
side-by-side at a site in Huntsville, AL. They estimated the sampling error for reflectivity 
(in dBZ units) as 1.36 dB (time window for SSD integration was 1-min). They also 
estimated the normalized sampling error for SR as 8.5%. The temporal correlation 
between the three measures of Z e is visually quite good. 
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By re-sampling the 2DVD and POSS reflectivities to the Kumpula radar samples, the 
scatter plot of 30 December 2010 case shown in Fig. 1 1 is obtained. The bias between 
POSS and 2DVD Z e is 0.11 dB, the standard deviation is 2.9 dB, and the correlation 
coefficient is 0.92. The corresponding values between Kumpula radar and 2DVD are, 
respectively, 0.18 dB (slight radar overestimate), 4.68 dB and 0.8. The latter standard 
deviation values would be even lower if the 2DVD sampling error of 1.3 dB were 
accounted for. 

The 2DVD processing described thus far gives the time series of Z e and SR every 
minute (i.e., 1-min time integration) for each of the long duration (> 4 h) events listed in 
Table 1. In order to realize a ‘stable’ Z e -SR relation the sequential intensity filtering 
technique (SIFT) described by Lee and Zawadzki (2005) is used along with weighted 
total least-squares to estimate the coefficient and exponent of the Z e -SR power law. The 
basic time window (W) selected is 1 h; the SSDs are ordered by increasing Z e in this 
window; and a moving average of M=5 consecutively ordered SSDs is done to filter the 
DSDs. The same procedure is performed for the next hour of the event and so on until 
the entire snow duration is covered. From the filtered DSDs, the Z e is re-computed using 
the appropriate p-D app power law. To re-compute SR, eq. (11) can no longer be used, 
rather it is computed as: 

SR = p(D app )Dl pp V f N{D app )dD app (13) 

where p = orD ap / is the mean fit, and Vf = c - d*exp(- KD app ) is the mean fit to the 
measured fall speeds (see example in Fig. 5). 
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Fig. 12 shows the scatter plot of Z e versus SR and the power law fit for (a) the entire 30 
Dec 2010 case (i.e., inclusive of both events listed in Table 3) and (b) 6 Jan. 2011 case. 
Table 3 shows the Z e -SR power law fits for the other three snow days. For reference the 
FMI climatological relation is Z e =100 SR 2 (Saltikoff et al. , 2010) It is fairly evident that for 
a given Z e , the SR from the FMI relation will exceed that predicted from Table 3 power 
law fits. For completeness Table 4 shows the ^-adjustment values, and the parameters 
[c d k] of the VrD a pp fit. 

3.3 Radar-derived snow accumulations 

There were three C-band polarimetric radars operating during LPVEx, being located at 
(see Fig. 4) Kumpula (KUM), Kerava (KER) and Vantaa (VAN). The technical 
specifications can be found in (Hirsikko et al, 2013; Saltikoff and Nevvonen, 2011). 
When radar reflectivity is used along with a Z e -SR relation to generate, e.g., daily (liquid 
equivalent) snow accumulation maps, clutter and beam-blockage at low elevation 
angles can cause loss of signal (in the case of clutter, due to clutter filtering) which 
manifests as artifacts in the snow accumulation maps. To avoid this problem, Z e data 
from the three radars have been composited, using maximum reflectivity factor from any 
of the three radars, to generate the snow accumulation map for 30 Dec 2010 as shown 
in Fig. 13. The peak accumulation is around 12 mm within the city of Helsinki. The solid 
black dots are the locations of six FMI snow gauges (Vaisalla VRG101 with Tretyakov 
wind shield; Lanza et al. 2006) and the OTT-PLUVIO at Jarvenpaa site. The numbers 
adjacent to the gauge locations are the measured accumulations in mm. The radar 
composite, of course, depicts quite clearly the spatial variability without any artifacts due 
to clutter or beam blockages; moreover the radar-based accumulations are in good 
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agreement with the gauges. Fig. 14 shows the accumulation map using the 
climatological FMI Z e -SR relation and it is readily apparent that, while the spatial 
variability is generally preserved, the magnitudes of the accumulations are 
overestimated relative to the gauges. In particular, the peak accumulations are now 
around 16 mm within the city. 


To further detail the radar and gauge comparisons, hourly accumulations are compared 
in Fig. 15 from one gauge location (solid dot in Fig. 13 with 7.8 mm; this is the Porvoo 
Harabacka location). Whilst it is clear that the FMI climatological relation overestimates 
the hourly accumulations soon after the snow begins, the radar-based hourly 
accumulation agrees well with the gauge (and not just the event totals). 


Fig. 16 (panel a) shows the scatter plot of daily gauge accumulations versus radar- 
based accumulations (extracted from the radar composites over the six gauge locations 
and the gauge at the Jarvenpaa site) for the 4 days using the Z e -SR power laws from 
Table 3 while panel b shows the same except for using the fixed FMI climatological Ze- 
SR relation. The significant feature is the dramatic reduction in bias resulting from using 
the Z e -SR obtained from 2DVD data as listed in Table 3 as compared with the fixed FMI 
relation. The normalized bias and normalized standard error values are, respectively, 
28% and 30.8% when Table 3 is used versus 96.6% and 66.1% for the fixed FMI 
relation. Note that positive bias implies radar overestimates the gauge values. The 
slope of a straight line trend passing through the origin is 1 .2 for panel a versus 1 .85 for 
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panel b. It is reasonable to infer that the FMI gauges could have underestimated the 
snow amounts due to wind and type of shielding (i.e., collection efficiency < 1). Recall 
that the FMI gauges are Vaisala VRG101 with Tretyakov wind shields whose collection 
efficiency is not known as a function of wind speed. The collection efficiency (or, 
undercatch) is a complicated function of not only gauge/shield type and wind speed, but 
also the type of snow particle (dry vs. wet or unrimed vs. rimed) and particle size 
distribution (Theriault et al. 2012). Thus, there is considerable scatter of the collection 
efficiency for a given wind speed along with a systematic decrease with increasing wind 
speed. The latter can be estimated from Rasmussen et al. (2012; their Fig. 1 1) as mean 
collection efficiency dropping to 0.75 at wind speed of 4 ms -1 . If this is taken into 
account the bias between radar and FMI gauges seen in Fig. 16a would be further 
reduced. 


3.4 Spatial reflectivity structure for 30 Dec 2010 and 06 January 2011 cases 

So far, the reflectivity structure nor the environmental/synoptic conditions have been 
described for the different snow days, as the main emphasis was on the 2DVD data, its 
processing and product evaluation. However, it is useful to consider the reflectivity 
structure for the 30 Dec 2010 case (which had the most daily accumulation) and the 06 
Jan 2011 case which had the least (Table 2 or Fig. 16), accompanied by very different 
coefficients/exponents of the Z e -SR power laws (respectively, 210/1.63 versus 
130/1.44). 
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On 30 Dec 2010 large scale snowfall areas from ESE (the first snowfall from 0800-1300 
UTC; see Fig. 10a) and from WNW (second snowfall from 1500-2400 UTC) merged 
above southern Finland. These snowfall areas were associated with two low pressure 
systems, one centered in Eastern Europe and the main one forming NW of 
Scandinavia. At around 1900 UTC the two precipitation systems have merged. 
ADMIRARI (Battaglia et al., 2010) LWP (liquid water path) observations reached a 
maximum of 400 g nr 2 at 1500 UTC. During the observations ADMIRARI was located in 
the backyard of Vaisala which is around 10 km north from Kumpula radar, as shown in 
Fig. 4. This is the time when the warm moist area from NW had arrived to the Helsinki 
region. During the snowfall the LWP values were ranging between 100-150 g rrr 2 . It 
should be noted that these are the slant LWP observations with elevation angle of 30°. 
At the time of peak snowfall for the first event (1100-1200 UTC; see Fig. 10a) SVI 
images were viewed and it was observed that the main precipitation types were pristine 
dendritic type crystals with large aggregates composed of dendrites (~ 8 mm) with little 
evidence of riming (Newman, personal communication). Fig. 17 shows sample SVI 
images at 1120 UTC near the time of maximum Z e (see Fig. 1 0a). 


On 06 Jan 2011 south westerly upper level flow from Scandinavia brought warm and 
moist air that resulted in a light to moderate snowfall lasting from 0100 to 0800 UTC. 
During this period (0500-0700 UTC) the ADMIRARI showed a large amount of 
supercooled water with LWP values exceeding 500 g rrr 2 . Examination of SVI images 
between 0600-0630 showed definite indications of rimed dendrites and columnar 
crystals followed by rimed to heavily rimed particles (perhaps graupel). Further in time, 
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large aggregates (~5-7 mm) appear to be rimed. Also, many smaller rimed 
snowflakes/crystals (Newman, personal communication). Fig. 18 shows examples of 
SVI images of rimed aggregates at 0620 UTC. 


Fig. 19 shows the rather dramatically different reflectivity structures (at low elevation 
angle 0.5°) between the 30 Dec 2010 event (at 1000 UTC) and the 06 January 2011 
event (at 0610 UTC). The spatial variability is much more pronounced in the 06 January 
case (more cellular) as compared to the more conventional spatial variability occurring 
on 30 Dec. The cellular feature implies weak imbedded convection is likely with more 
prevalent particle riming as alluded to earlier. This is supported by analysis by Lim et al. 
(2013) who were able to associate higher spatial variability with enhanced riming of 
particles. The vertical structures are also different as depicted in the RHI scans in Fig. 
20 taken along the radial to the Jarvenpaa site Again, the 30 Dec event is much more 
uniform in the vertical as compared with the 06 Jan event, the latter showing more 
evidence of cellular structure in the vertical implying imbedded convection and 
enhanced riming. 


Finally, in Fig. 21, the snow accumulation map for 06 Jan event is shown using the 
2DVD-derived Z e -SR power law, which can be compared with Fig. 13. The snow 
accumulations for this event (see, also Table 1) are much smaller compared with 30 
Dec case additionally showing much more spatial variability. Note these are not daily 
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totals but restricted to the period 0000-0900 UTC since the 2DVD stopped working at 
0824 UTC on this day thereby missing another major snowfall event later on this day. 

4. Summary and Conclusions 

The estimation of the mean density-size and Z e -SR power laws using 2D-video 
disdrometer measured fall speed, apparent diameter and snow size distribution (SSD) 
along with Bohm’s (1989) methodology is described in some detail. A method for 
adjusting the concentration based on single camera data to account for loss of particles 
that do not satisfy the matching criteria (when 2 cameras are used) is shown to be 
reasonable when compared with Snow Video Imager (SVI)-based concentrations. 
Snow events which occurred on four days of the Light Precipitation Validation 
Experiment (LPVEx) were chosen based on light wind speeds (< 4 ms -1 ) at the 
measurement site with liquid equivalent snow accumulations ranging from 1 .5 mm to 4 
mm. While there is large variability of fall speed, area ratio and derived density which is 
attributed to natural variability of snow type, shape and porosity, the mean density-D app 
(or, mass-D apP ) and Z e -SR power laws do vary from event-to-event. The reflectivity 
derived from the 2DVD data is shown to be in good agreement with collocated POSS 
and with scanning C-band radar, while the liquid equivalent snow accumulation is 
shown to be in good agreement with collocated OTT-PLUVIO gauge at the 
measurement site. The radar-based snow accumulations using the 2DVD-derived Ze- 
SR relations for the four days are in good agreement with a network of six FMI snow 
gauges (and the OTT-PLUVIO gauge) and outperform the accumulations derived from a 
climatological Z e -SR relation used by the Finnish Meteorological Institute (FMI). The 
normalized bias between radar-derived and gauge accumulation is reduced from 96% 
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(overestimate by the FMI climatological relation) to 28% when using the Z e -SR based 
on 2DVD data. The normalized standard error is also reduced significantly from 66% to 
31%. While the FMI gauges were equipped with Tretyakov wind shields, 
undercatchment due to wind cannot be ignored and could account for underestimation 
of snow accumulations by 20-30% for wind speeds in the range 3-4 ms -1 ; this would 
reduce the bias between radar and gauge accumulations even further. 


For two of the days with widely different coefficients of the Z e -SR power law, the 
reflectivity structure showed significant differences in spatial variability (both horizontal 
and vertical). Liquid water path estimates from radiometric data also showed significant 
differences between the two cases. Examination of SVI particle images at the 
measurement site corroborated these differences in terms of unrimed versus rimed 
snow particles. 


In summary, the findings reported herein support the application of Bohm’s (1989) 
methodology for deriving the mean density-size and Z e -SR power laws using data from 
2D-video disdrometer. Evaluation of radar-based snow accumulation against a network 
of snow gauges independently supports the latter conclusion notwithstanding the limited 
number of events available for analysis. 
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List of Figures 


Fig 1 : Front- and side-view images of a snow flake from the 2DVD. The maximum ‘size’ 
is approximately 16 mm in the side view (Camera B); from Huang et al. (2010). 


Fig. 2: Illustrating the measurement principle of the 2D-Video disdrometer (from Kruger 
and Krajewski 2002; Schonhuber et al. 2008). Note that the virtual measuring area is 
different (100 cm 2 ) relative to the single camera measuring area (250 cm 2 ). 


Fig. 3: An example to show the shadow area ( A e ) and the smallest circumscribed area 
area (A); from Bohm (1989). 


FIG. 4: The map shows the location of 3 C-Band radars (KUM, VAN and KER), the 
Jarvenpaa site and other 6 FMI gauges used in this paper. Note that the OTT-PLUVIO 
gauge and 2DVD are sited at Javenpaa. 


FIG. 5: Example of 2DVD-derived fall speeds versus apparent diameter for snow event 
from 30 Dec 2010 dominated by snow aggregates. For each size bin the mean and ±1o 
standard deviation of the fall speed are shown. The exponential fit to the mean values of 
the form Vf = 1 .12-1 ,31*exp(-1 .42 D app ) where D app is in mm and Vf in m s _1 is also 


shown. 
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Fig. 6: Particle (snow) size distributions from the 2DVD and from the snow video imager 
(SVI) for (a) 30 Dec 2010 first snow event with y equal to 2.21 (b) 6 Jan. 2011 event 
with y equal to 2.55. 


Fig. 7: Frequency of occurrence (in color log scale) of A- vs. D app for 30 Dec 2010 snow 
event; also shown are the bin averaged mean and ±1a standard deviation bars along 
with the power law A,-D app fit. 


Fig. 8: Density versus D app (mean and ±1a) along with mean p-D app power law fit for the 
30 Dec 2010 event. 


Fig. 9: Liquid equivalent snow accumulation for 30 Dec 2010 comparing 2DVD-derived 
accumulations against collocated OTT-PLUVI02 gauge. 


Fig. 10: Time series comparison of reflectivity from 2DVD, POSS and scanning 
Kumpula radar for (a) 30 Dec 2010 event and (b) 6 Jan. 2011. 
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Fig. 11: Scatter plot of 2DVD or POSS reflectivity (Y-axis) versus Kumpula radar 
reflectivity (X-axis) extracted over Jarvenpaa site for 30 Dec 2010 event. 


Fig. 12: Z e versus SR scatterplot using SIFT method for (a) 30 Dec 2010 case and (b) 6 
Jan. 201 1 . The power law fit using weighted total least squares is also shown. 


Fig. 13: Snow accumulation map for 30 Dec 2010 (entire day) using Z e -SR power law in 
Table 3. Reflectivity data from three radars have been composited. The solid black dots 
are locations of 6 FMI snow gauges and the numbers adjacent are the snow gauge 
measurements of liquid equivalent snow accumulation in mm. 


Fig. 14: As in Fig. 13 except the FMI climatological Z e =100*SR 2 is used. 


Fig. 15: Comparison of hourly accumulations at the Porvoo Harabacka site using FMI 
gauge. 


Fig. 16: Scatter plots of daily gauge accumulations from the seven FMI snow gauges 
versus radar-based daily accumulations using (top panel a) daily Z e -SR relations from 
Table 3 and (bottom panel b) using fixed FMI climatological Z e -SR relation. Data from 
the 4 snow days are color-coded. 
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Fig. 17: Example of SVI images of large aggregates at 1120 UTC on 30 Dec 2010 at 
Jarvenpaa site (also close to time of peak reflectivity, see Fig. 10a). 


Fig. 18: Example of SVI images of large rimed aggregates at 0620 UTC on 06 Jan 201 1 
at Jarvenpaa site. 


Fig. 19: Spatial variability of reflectivity (left panel) for 30 Dec 2010 and, (right panel) for 
06 Jan cases. 


Fig. 20: RHI scan data from Kumpula radar along the radial to the Jarvenpaa site (top 
panel) 30 Dec 2010 at 1053 UTC, and (bottom panel) 06 Jan 201 1 at 0603 UTC. 


Fig. 21 : As in Fig. 13 except for 06 Jan 201 1 event (0000-0900 UTC only) using 2DVD- 


derived Z e -SR relation. 
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Fig. 1: Front- and side-view images of a snow flake from the 2DVD. The maximum 
‘size’ is approximately 16 mm in the side view (Camera B); from Huang et al. (2010). 
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Fig. 2: Illustrating the measurement principle of the 2D-Video disdrometer (from Kruger 
and Krajewski 2002; Schonhuber et al. 2008). Note that the virtual measuring area is 
different (100 cm 2 ) relative to the single camera measuring area (250 cm 2 ). 
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Fig. 3: An example to show the shadow area (A e ) and the smallest circumscribed area 
area (A); from Bohm (1989). 
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Fig. 4: The map shows the location of 3 C-Band radars (KUM, VAN and KER), the 
Jarvenpaa site and other 6 FMI gauges used in this paper. Note that the OTT-PLUVIO 
gauge and 2DVD are sited at Javenpaa. 
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30 December 2010; 00:00-15:00 UTC 



app 

Fig. 5: Example of 2DVD-derived fall speeds versus apparent diameter for snow event 
from 30 Dec 2010 dominated by snow aggregates (00:00-1 5:00 UTC). For each size bin 
the mean and ±1o standard deviation of the fall speed are shown. The exponential fit to 
the mean values of the form Vf = 1 .12-1 ,31*exp(-1 .42 D app ) where D app is in mm and Vf 
in m s' 1 is also shown. 


30 December 2010; 00:00-15:00 (UTC) 
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Fig. 6: Particle (snow) size distributions from the 2DVD and from the snow video imager 
(SVI) for (a) 30 Dec 2010 first snow event with y equal to 2.21 (b) 6 Jan. 2011 event 
with y equal to 2.55. 
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Fig. 7 ; Frequency of occurrence (in log scale) of A r vs. D app for 30 Dec 2010 snow 
event; also shown are the bin averaged mean and ±1 a standard deviation bars along 
with the power law Ar-D app fit. 
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Fig. 8: Density versus D app (mean and ±1 a) along with mean p-D app power law fit for the 
30 Dec 2010 event. 
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Fig. 9: Liquid equivalent snow accumulation for 30 Dec 2010 comparing 2DVD-derived 
accumulations against collocated 0TT-PLUVI02 gauge. 
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Fig. 10: Time series comparison of reflectivity from 2 DVD, POSS and scanning 
Kumpula radar for (a) 30 Dec 2010 event and (b) 6 Jan. 2011. 


30 



Fig. 11: Scatter plot of 2DVD or POSS reflectivity (Y-axis) versus Kumpula radar 
reflectivity (X-axis) extracted over Jarvenpaa site for 30 Dec 2010 event. 
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Fig. 12: Z e versus SR scatterplot using SIFT method for (a) 30 Dec 2010 case and (b) 6 
Jan. 2011. The power law fit using weighted total least squares is also shown. 
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Fig. 13: Snow accumulation map for 30 Dec 2010 (entire day) using Z e -SR power law in 
Table 3. Reflectivity data from three radars have been composited. The solid black dots 
are locations of 6 FMI snow gauges and the numbers adjacent are the snow gauge 
measurements of liquid equivalent snow accumulation in mm. 
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FIG. 14: As in Fig. 13 except the FMI climatological Z e =100*SR 2 is used. 


Hourly precip accumulation - Porvoo Harabacka (WMO 02983) 



Fig. 15: Comparison of hourly accumulations at the Porvoo Harabacka site using FMI 
gauge. 
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Fig. 16: Scatter plots of daily gauge accumulations from the seven FMI snow gauges 
versus radar-based daily accumulations using (top panel a) daily Z e -SR relations from 
Table 3 and (bottom panel b) using fixed FMI climatological Z e -SR relation. Data from 
the 4 snow days are color-coded. 


834 




3463 fpm 2010 12 30 11 20 p. 1 


5 mm 



835 


Fig. 17: Example of SVI images of large aggregates at 1120 UTC on 30 Dec 2010 at 
Jarvenpaa site (also close to time of peak reflectivity, see Fig. 10a). 
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Fig. 18: Example of SVI images of large rimed aggregates at 0620 UTC on 06 Jan 
2011 at Jarvenpaa site. 
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Fig. 19: Spatial variability of reflectivity (left panel) for 30 Dec 2010 and, (right panel) for 
06 Jan cases. 
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Fig. 20: RHI scan data from Kumpula radar along the radial to the Jarvenpaa site (top 
panel) 30 Dec 2010 at 1053 UTC, and (bottom panel) 06 Jan 2011 at 0603 UTC. 
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Fig. 21: /As in Fig. 13 except for 06 Jan 2011 event (0000-0900 UTC only) using 2DVD- 
derived Z e -SR relation. 
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852 Table 1 : The four snow days from LPVEx 


Event 

Time (UTC) 

Temperature 

C 

2DVD 

accum 

(mm) 

0TT-PLUVI02 accum 
(mm) 

30 Dec 2010 

a. 0800-1 300 

b. 1500-2359 

a. -9 

b. -8 

3.84 

4.24 

06 Jan 2011 

0230-0830 

-7 

1.53 

1.53 

12 Jan 2011 

a. 0800-1230 

b. 2230-2359 

a. -3 

b. -3 

3.36 

2.05 

13 Jan 2011 

0000-0500 

-4 

2.77 

2.73 


853 


Table 2: Coefficient o and exponent (/ 3 ) of p-D app power law fit (density in g cnr 3 and 
Dapp in mm) 


Event 

Time (UTC) 

a 

P 

30 Dec 2010 

a. 0800-1 300 

b. 1500-2359 

a. 0.15 

b. 0.15 

-0.86 

-0.96 

06 Jan 2011 

0230-0830+ 

0.17 

-0.79 

12 Jan 2011 

a. 0800-1230 

b. 2230-2359 

a. 0.23 

b. 0.19 

-0.88 

-0.8 

13 Jan 2011 

0000-0500 

0.19 

-0.8 


f the 2DVD stopped working at 0824 on this day 
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Table 3: The Z e = a*SR b power law for the four days. Note Z e in mm 6 nr 3 and SR in mm 
h' 1 . 



a 

b 

30 Dec. 2010 

210.72 

1.63 

06 Jan. 2011 

130.72 

1.44 

12 Jan. 2011 

209.20 

1.67 

13 Jan. 2011 

134.86 

1.81 


855 f The power laws are derived for application to all events occurring during the 

856 day. 

857 


858 


Table 4: The ^adjustment factor and the Vf-D app fit parameters [c d k]. 


Event 

Time (UTC) 

r 

ysvi 

c 

d 

K 

30 Dec 2010 

a. 0800-1 300 

b. 1500-2359 

a. 2.21 

b. 1.94 

a. 2.46 

b. 2.57 

1.20 

1.54 

1.16 

06 Jan 2011 

0230-0830 

2.55 

2.92 

1.37 

1.37 

1.25 

12 Jan 2011 

a. 0800-1230 

b. 2230-2359 

a. 3.53 

b. 2.23 

a. 3.6 

b. 2.18 

1.37 

1.85 

1.78 

13 Jan 2011 

0000-0500 

2.23 

4.24 

1.34 

1.34 

0.95 


f The Vf- Dapp fits are derived for application to all events occurring during the day 
(note: D app in mm and Vf in m s' 1 ). 


